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In the corner-sharing lattice, magnetic frustration causes macroscopic degeneracy in the 
ground state, which prevents systems from ordering. However, if the ensemble of the degenerate 
configuration has some global structure, the system can have a symmetry breaking phenomenon 
and thus posses a finite temperature phase transition. As a typical example of such cases, the 
magnetic phase transition of the Ising-like Heisenberg antiferromagnetic model on the kagome 
lattice has been studied. There, a phase transition of the two-dimensional ferromagnetic Ising 
universality class occurs accompanying with the uniform spontaneous magnetization. Because 
of the macroscopic degeneracy in the ordered phase, the system is found to show an entropy- 
driven ordering process, which is quantitatively characterized by the number of "weathervane 
loop" ( WL) . We investigate this novel type of slow relaxation in regularly frustrated system. 
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In the past several decades, there has been increas- 
ing interest in spin structure of several frustrated spin 
system, where non-coUinear spin configurations appear 
due to competition of interaction. The structure and 
symmetry of the non-coUinear structures cause various 
types of interesting ordering processes.^' ^ In particu- 
lar, antiferromagnets on the so-called corner-sharing lat- 
tices, e.g., kagome and pyrochlore lattices, the frustra- 
tion causes macroscopic degeneracy even in continuous 
spin systems, and thus in contrast to the case of trian- 
gular lattice, magnetic ordering in XY and Heisenberg 
models is difficult in the kagome lattice. However, for 
classical Heisenberg antiferromagnet on the kagome lat- 
tice, Harris et al.^ estimated magnetic susceptibility by 
means of high temperature expansion up to ffi , and con- 
cluded that -x/S X -s/S structure should be selected for 
XY and Heisenberg spin system. Huse and Rutenberg"' 
and Chalker et al.^ studied hidden order parameter at 
T 0+. Moreover, the concept of order by disorder 
has been discussed to choose an ordered state in en- 
ergetically degenerate states. Experimentally, various 
kagome antiferromagnetic compounds have been stud- 
ied, e.g., jarosite-based materials AM3(OH)6(S04)2 (A+ 
is a caution such as K+, Rb"*", NH4 , T1+, Ag+ and 
Na+, etc. and is a magnetic ion such as Fe'^+ or 



Cr3+ ete.),i"-i5 SrCr9,Gai2_9,Oi9,i6.i7 Rb2M3S4 {M^+ 
is a magnetic ion such as Ni^+ , Co^+ , Mn^+ ),^^ 
Zn^Cu4-:,(OH)6Cl2,^^ Ba2Sn2Ga3ZnCr7022^° and ^Hc 
on the graphite sheet^-'^ etc. 

Beside the isotropic Heisenberg model, Kuroda and 
Miyashita^^ pointed out that a thermodynamical phase 
transition takes place in the Ising-like Heisenberg kagome 
antiferromagnets. The properties of this phase transition 
were also recently investigated in details by Bekhechi 
and Southern.^"* There, the ground state magnetization 
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has two-fold degeneracy, and a phase transition which 
breaks this symmetry occurs at finite temperature. Be- 
cause of the nature of antiferromagnets on corner sharing 
lattices, there still remain macroscopic degenerate states 
in the magnetically ordered state. This macroscopic de- 
generacy causes peculiar properties in the ordered state 
such as the localization of the spin wave, etc. and the 
low temperature magnetic state is called "exotic ferro- 
magnetic phase". As to experimental realization of this 
phase transition, Maegawa et al. found two cusps of the 
temperature dependence of the susceptibility in kagome 
system NH4Fe3(OH)6(S04)2,^^ and they concluded that 
the successive phase transition may be caused by an easy 
axis type anisotropy in the Heisenberg kagome antiferro- 
magnet. 

The slow relaxation has been studied in the spin- 
glasses and also diluted ferromagnets, where the random- 
ness of the system plays an important role for the slow 
dynamics. Recently, however, slow relaxation processes 
have been observed also in regularly frustrated spin sys- 
tems such as (H30)Fe3(OH)e (804)2^^ experimentally. 
Generally we may expect that the frustration prevents 
ordering and should cause the fast dynamics. However, 
we pointed out that there exists slow relaxation in an 
entropy-driven ordering process in macroscopically de- 
generate states (i.e., a kind of the order by disorder mech- 
anism) in the Ising-like Heisenberg kagome antiferromag- 
netic model. Bekhechi and Southern^"' also found slow 
dynamics and discussed it in the context of the glassy dy- 
namics. In this Letter, we investigate microscopic proper- 
ties of this slow relaxation by studying structure change 
in macroscopically degenerate states which are charac- 
terized by the "weathervane loop". We find an Arrhe- 
nius type relaxation process in the magnetically ordered 
state. The relaxation time of this slow process is much 
longer than that of magnetization. 

We consider Ising-like Heisenberg antiferromagnetic 
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model in the kagomc lattice. The Hamiltonian is given 

by 

7^ = - J ^ {SfSJ + SfS^ + AS^S^) , (A > 1) , (1) 

where {i, j) denotes nearest neighbor pairs in the kagome 
lattice and A denotes Ising-like anisotropy. 

The kagome lattice has a corner-sharing struc- 
ture of triangle clusters and the spins on a tri- 
angle cluster consists of the set {Sa,Sf3,S-y) in the 
ground state. For the Ising-like Heisenberg anti- 
ferromagnetic model, these are given by Sa ~ 
(0,0,1), Sp — (sin 6* cos 0, sin 6* sin 0, — cos 6*) and ~ 
(— sin0cos0, — sin^sin^, — cos9), where cos9 = A/(l-f 
A). The angle (j) is any angle denoting the freedom of ro- 
tation of 2tt. 

There are two types of ground state configurations at 
each triangle. One is the one mentioned above, which we 
call "up-down-down" configuration. The other is the one 
of up-side-down structure, that is Sa = (0,0,-1), etc., 
which we call "down-up-up" configuration. Although the 
positions of the spins {Sa,Sp,S-y) is not fixed, the ground 
state has two- fold degeneracy. That is, if one triagle has 
the "up-down-down" configuration, all others must have 
the same type. Because "up-down-down" configuration 
has a nonzero negative magnetization, and the "down- 
up-up" configuration has the opposite one, this symme- 
try breaking causes appearance of the spontaneous mag- 
netization. That is, all the triangles have the local mag- 
netization M = ±(1 — 2 cos 6), and the symmetry break- 
ing of the magnetization takes place, which causes the 
thermodynamic phase transition at finite temperature 
where Z2 symmetry is broken in spite of no sublattice 
long range order exists. ^'^ 

If we connect the spins of Sf} or on the lattice, 
closed loops arc formed as depicted by the bold red lines 
in Fig. 1. These closed loops are called "weathervane 
loop"(WL). There are macroscopic ways to draw the 
WLs on the lattice. Indeed the number of the ways is 
given by Udagawa et al. which is the same as that of the 
S = 1/2 Ising antiferromagnetic kagome model in the 
magnetic field. ^® All of them are energetically degener- 
ate. 



rotation 27r at each WL, which causes the residual en- 
tropy. We investigate dynamics of the number of the WLs 
'T'loop to characterize the relaxation of spin structure. 
The maximum number of WLs N^^^, is realized in the 

\/3x^/3 structure, in which all WLs are a simple hexagon 
which has minimum length. In this case, N^^^ is 
where N is the number of sites. At finite temperature, 
there exist broken WLs with energetic defects. We call 
these broken loops "weathervane strings". We studied 
relaxation processes from a random configuration by the 
thermal bath Monte Carlo method in the lattice with 
N = 675. In the present work, we adopted A = 3. We 
obtained the data by averaging over 24 independent runs. 
We depict the dynamics of these quantities at T = 0.05 J 
which is below (= 0.076 J ± O.OOIJ)^'* in Fig. 2. 

When the energy or the magnetization (red points in 
Fig. 2) reaches to its equilibrium value, the total length 
of the WLs and strings (blue points in Fig. 2) also reaches 
to the equilibrium value. Here the magnetization means 
the average of the z-componcnt of spins. However, the 
ratio of the contribution of the loops (green points in 
Fig. 2) increases in time as depicted in the Fig. 2. 
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Fig. 2. (Color online) Relaxation of nioop, string and magnetiza- 
tion at T = 0.05 J from random configuration. 



Figures 3 are the snapshots of spin configuration at 
10^, 10"^, and 4 X 10"^ MCS. The red closed and green open 
lines denote WLs and weathervane strings, respectively. 




Fig. 1. (Color online) Spins on each triangle cluster consist of the 
three types {Sa, S/j, S-y). The thick line denotes the WL. 



In macroscopic energetically degenerate states the ef- 
fect of spin configuration entropy sometimes plays an im- 
portant role to determine the thermodynamical ordered 
state. In the present model, there exists the freedom of 




Fig. 3. (Color online) Snapshots at lO^MCS, lO'^MCS, 4 X 
lO^MCS. The bold lines denote WLs and strings, respectively. 
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Let US study temperature dependence of the entropy of 
a WL of the length L at finite temperatures. At low tem- 
peratures, the low energy excitation is restricted in the 
freedom of the angle of {Sp} and {S^}. This degree of 
freedom is expressed by the antiferromagnetic XY chain 
with periodic boundary condition. Ti. = —J^i^ij) • Sj, 
where Si = (cos0, sin0). The entropy per spin decreases 
with the length and the effect of spin chain length be- 
comes more significant when the temperature decreases. 
Because the entropy per spin decreases monotonously 
with the length of WL, the short loops are entropically 
preferable, and thus nioop should increase at low temper- 
atures. What is more, at T ^ 0+, we expect that nioop 
must be the maximum value and the spin structure be- 
comes the so-called -y/S x -y/S structure. 

Next, we study the relaxation of magnetization and 
"■loop- We ready the three types of initial configurations, 
i.e. (a) the \/3 x ^/i structure, (b) the q — Q structure 
and (c) a random structure. The configurations (a) and 
(b) are typical ground states of the present model, and 
the configuration (c) corresponds to a state just after 
quench the temperature from a high temperature. 
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case (a) at various temperatures below the critical tem- 
perature. There, we find two steps in relaxation. The first 
relaxation corresponds to initial local relaxation from the 
complete ground state configuration. The relaxation time 
of this process is very short and does not strongly de- 
pend on the temperature. In the second relaxation, the 
re-construction of the WLs takes place. The first and sec- 
ond relaxations correspond to energetical and cntropical 
relaxations, respectively. It is noted that the two-step 
relaxation is a characteristics of the ordered state of the 
present model. 
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Fig. 5. (Color online) Relaxation of nioop from \/3x structure 
at several temperatures. Dashed lines denote the fittling curves 
estimated by Eq. (2). 



Fig. 4. (Color online) Relaxation of the magnetization and nioop 
at T = 0.05J from (a) \/3 X configuration, (h) q = config- 
uration and (c) random configuration. 

In Figs. 4, the relaxation processes at T = 0.05 J are 
plotted. In the cases (a) and (b), the magnetization is 
maximum at t = 0, and it relaxes very fast to uniformly 
magnetized ordered state. The relaxation of magnetiza- 
tion to the equilibrium is depicted in the inset. In con- 
trast, in the case (c), i.e., from a random state, it takes 
some time to realize the uniformly magnetized state. 
Thus we regard the relaxation time in the case (c) as 
the intrinsic relaxation time of the magnetization Tmag- 

The relaxation of nioop is much slower. In the case (a) 
nioop starts with the maximum value, while in the case 
(c) it starts with zero. In all the cases, (a), (b) and (c), the 
numbers seem to reach the same saturated value. This 
observation indicates that there exists stable equilibrium 
state for nioop ■ The relaxation time rioop is also about the 
same, although in the case (b) a non-monotonic process is 
observed in the early state where the weathervane lines 
in the initial q = Q state are broken into short strips. 
Thus, we expect that there exists an intrinsic relaxation 
time for nioop- From the figure we estimate that Tmag is 
about lO'' MCS, and rioop is of the order lO^MCS. 

In Fig. 5, we plot the relaxation process of nioop in the 



We fit all the second relaxation processes in the fol- 
lowing form 

nioop {t) = A (T) exp (-7^) + «S (T) , (2) 

T (T) denotes a characteristic relaxation time of the 
'^loop and n|^^p is the equilibrium value of nioop- Because 
the first relaxation is much faster than the second one, 
the choise of start time of the second relaxation is irrel- 
evant. Here we fit the curves in the cases of T = 0.07 J, 
0.065J, 0.06J, 0.055J, 0.05J, 0.0475J, 0.045J, 0.0425J 
and 0.04J. The temperature dependence of r and nj^^p 
are plotted in Figs. 6. In Fig. 6(a), we find that nioop in- 
creases when the temperature is lowered. We expect that 
this value continues to the ground state value nioop = 1- 
It should be noted that this quantity is not zero at the 
critical temperature, and it can not be an order param- 
eter although it is a good indicator of the degree of the 
order of the structure. In Fig. 6(b), we find a good lin- 
earity, which indicates the Arrhenius law t(T) cx e^^^ . 
From the slope in Fig. 6(b), the energy barrier /S.E is 
estimated about 0.78 J ~ O (J). It is expected that this 
energy is necessary when the WLs are reconnected lo- 
cally. 

We study slow relaxation of spin configuration in the 
magnetically ordered phase of the Ising-like Heisenberg 
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Fig. 6. (Color online) Temperature dependencies of equilibrium 
values of magnetization (closed circles for A'^ = 675, open circles 
for A'^ = 972, closed diamonds for N = 1728, open diamonds for 
N = 2700, respectively). Temperature dependence of equilibrium 
values of ni^op estimated by Eq. (2) (crosses). The inset shows 
the relaxation time t as a function of inverse temperature /3. The 
dashed line denotes r = O.SSe"'^*''. 



kagome antiferromagnets. From the view point of the 
entropy, states with the larger nioop are preferable, al- 
though macroscopic configurations are energetically de- 
generate and they have the same magnetization. It is ex- 
pected that the spin structure becomes -v/Sx structure 
atT ^0+ hmit. 

We study the relaxation of the number of the weath- 
ervane loops which characterizes the spin configuration, 
and we found the relaxation obeys the Arrhenius law r oc 
g/^AB ^[^Ij characteristic energy barrier AE ^ O (J). 

In the present work, we studied relaxation processes 
in the system of iV = 675, and pointed out two dif- 
ferent relaxation processes exist. They are energetical 
and entropical relaxations, respectively. However, we do 
not studied the size-dependence of the relaxation time. 
The energetical relaxation is a standard one of the fer- 
romagnetic model where we expect the y/i scaling, i.e. 
Tmag ~ iV (= L^) . On the other hand, the latter relax- 
ation time depends on how the reconnection of WLs take 
place, and at this moment we do not have definited con- 
clusion yet, which will be studied in near future. 

Finally, we would like notify that the present phenom- 
ena are the effect of residual entropy, but not the en- 
tropy of the excited levels which plays important roles 
in the reentrant phase transition. It is cither the en- 
tropy effect due to the simple harmonic potential which 
often plays important role to determine the ground 
state configuration from continuously degenerate ground 
state.30'31 
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